The Nature of Asymmetry in Fluid Criticality 
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By combining accurate liquid-vapor coexistence and heat-capacity data, we have unambiguously 
separated two non-analytical contributions of liquid-gas asymmetry in fluid criticality and showed 
the validity of "complete scaling" [Fisher et al, Phys. Rev. Lett. 85, 696 (2000); Phys. Rev. 
E, 67, 061506, (2003)]. We have also developed a method to obtain two scaling-field coefficients, 
responsible for the two sources of the asymmetry, from mean-field equations of state. Since the 
asymmetry effects are completely determined by Ising critical exponents, there is no practical need 
for a special renormalization-group theoretical treatment of asymmetric fluid criticality. 
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A controversial issue of liquid-gas asymmetry in fluids 
has been a subject of prolonged discussions for more than 
a century since Cailletet and Mathias' discovery of the 
empirical "law" of rectilinear diameter yj. According 
to this " law" , the mean of the densities of liquid p' and 
saturated vapor p" is a linear function of the temperature 
T: 



^Pd = -^ = 1 



D 



AT 



(1) 



where p^ is the critical density, AT = (T — Tc)/Tc the 
reduced distance to the critical temperature Tc- The 
system-dependent coefficient D generally increases from 
0.02 for •^He 2] to values larger than unity with increase 
of Tc 3, 4] . In this Letter we demonstrate how the mean- 
field rectilinear diameter splits up in the critical region 
into two " singular diameters" associated with two differ- 
ent sources of asymmetry. It is commonly accepted that 
fiuids asymptotically (AT -^ 0) belong to the critical- 
point universality class of the Ising model 5] . We argue 
that the asymmetry effects in near-critical fluids, at least 
in lower approximation, are determined by Ising criti- 
cal exponents, hence, in contrast to a commonly used 
approach, there is no need for a special renormalization- 
group theoretical treatment of asymmetric fluid critical- 
ity. 

Thermodynamics near a critical point is controlled by 
two scaling fields, "ordering" hi and "thermal" ft.2, while 
an appropriate field-dependent potential /13 is a universal 
function oi hi and /12 |5|: 



hi-v' 



hi 



(2) 



where a = 0.109, /3 = 0.326, and 7 = 1.239, interrelated 
as a -1-2/3 -|- 7 = 2, are universal Ising critical exponents in 
the scaling power laws (as a function of /12 at hi = 0) for 
the weakly-divergent susceptibility, order parameter, and 
strongly-divergent susceptibility, respectively |5|. The 



scaling " densities" , a strongly fluctuating order parame- 
ter (l)i and a weakly fluctuating (j)2, conjugate to hi and 
/i2, such that dh^ — c/fidhi + (j)2dh2. The universal func- 
tion / contains two system-dependent amplitudes, and 
the superscript ± refers to /12 ^ 0. The Ising model for- 
mulated for fluids is known as "lattice-gas" |g]. In the 
lattice-gas /13 is the "critical part" of the grand ther- 
modynamic potential, O = —PV, taken per unit volume, 
thus hi is the dimensionless chemical-potential difference 
Ap — (p — Pc)/k^Tc, where p^ is the value of the chemi- 
cal potential at the critical point and fee is Boltzmann's 
constant. By definition, hi — Q along the critical iso- 
chore above T^ and along the liquid-vapor coexistence 
curve below Tc. Hence, in the lattice gas the order pa- 
rameter (f)^ = Ap = (p — p^)/p^, the thermal scaling field 
/i2 = AT, and 02 = A{pS) = (p5 - pJc)/kB, where pS 
is the density of entropy. 

The lattice gas has a special symmetry: the order pa- 
rameter is symmetric with respect to the critical isochore 
{D = 0). Since the early 1970's the fiquid-gas asym- 
metry has been commonly incorporated into the lattice- 
gas analogy by linear mixing of two independent physical 
fields Ap, and AT into the both theoretical scaling fields 
hi and /12 .7|- Since the absolute value of entropy is 
arbitrary, mixing of AT into hi plays no role. Contrar- 
ily, mixing of Ap (with 62 as a mixing coefficient) into 
/i2 has an important consequence, known as the "sin- 
gular diameter" : the mean of the densities must con- 



tain a non-analytic contribution cc &2 



AT 



so that 



d{Ap^)/dT oc 62 AT , diverging weakly. However, 

the chemical potential would remain an analytical func- 
tion of temperature along the liquid-vapor coexistence 
and the symmetry would be restored by a redefinition of 
the order parameter as (l)i = Ap + b2A[pS). 

At this point we encounter a major problem. First of 



all, the existence of the 



AT 



term in the "diameter" 



of real fluids has never been detected unambiguously . 
While some fluids show strong deviations from rectilin- 
ear diameter, apparently even stronger than AT 



many fluids show very little or no deviations at all |a, |9j . 
Moreover, there is a conceptual problem with mapping 
real fluids into the lattice-gas even at the mean-fleld level. 
In the mean-field approximation the critical part of the 
thermodynamic potential is represented by a Landau ex- 
pansion: 



hs = -^00^201 + tMo0i - hi(j)i- 



(3) 



When hi = A/i, ft.2 = AT +b2Afi, and cj)^ = Ap + 
&2A(pS'), this expansion generates asymmetric terms oc 
b2AT{Ap)^ and oc b2{Ap)^. However, in the simplest 
equation of state that describes real-fluid behavior, the 
van der Waals equation, the term a AT{Ap)^ is absent, 
while the term oc (Ap)^ exists. Furthermore, in most 
classical equations of state dp? /dT^ along the liquid- 
vapor coexistence exhibits a discontinuity directly related 
to the existence of the independent 5th-order term [l^ 
in Landau expansion. The existence of the independent 
5th-order term makes exact mapping of fluids into the 
lattice-gas model by the conventional mixing of physical 
fields impossible. On the other hand, a renormalization- 
group treatment of the 5th-order term resulted in an in- 
dependent critical exponent 0^ ~ 1-3 JjJ (not existing in 
the Ising model!). We show, however, that asymmetric 
fiuids can be consistently mapped into Ising criticality 
by applying so-called "complete scaling" originally pro- 
posed by Fisher et al. [ij]. A redefinition of the order 
parameter, suggested by complete scaling, results in elim- 
ination of the 5th-order term in Landau expansion, thus 
making the renormalization-group treatment of the 5th- 
order term irrelevant for fluids. 

Complete scaling suggests that all three physical fields 
A/i, AT, and AP = (P - Pc)/Pc^bTc are equally mixed 
into three scaling fields /ii, ft.2, and ^13. In linear approx- 
imation 



hi 


= aiAjj, + a2AT + a^AP , 


(4) 


h2 


^ biAf + b2Ap + b3AP, 


(5) 


hz 


= CiAP + CsA/i + CgAf. 


(6) 



Before we apply complete scaling to describe asymme- 
try in fluids, we note that the relations between scaling 
and physical flelds can be simplified. The coefficients 
ai and 61 can be absorbed in two amplitudes in the 
scaling function /*. The coefficients ci and C2 are ab- 
sorbed by making the thermodynamic potential h^ di- 
mensional. The coefficient C3 = Sc is determined by 
the choice of the critical value of entropy. By adopt- 

ing ^c = (A^BPe)"' (dP/dT)^^^, ^ = (dP/df) , the 

\ / cxc,c 

slope of the saturation-pressure curve at the critical 
point, one obtains 02 = —03 (dP/dTj and 63 — 0. 

V / cxc.c 



Finally, the scaling fields contain only two amplitudes 
responsible for asymmetry in fiuid criticality: 



hi = Ap + asiAP- [dP/dT 

h2 = Af + b2Ap, 

hs = 



-AP + Ap+ (dP/df 



^J. J, 

cxc.c 


\' 1 




(8) 


AT. 


(9) 



As a result, while the order parameter in fiuids is, in 
general, a nonlinear combination of density and entropy 
01 = [Ap + b2A{pS)\/[l + asAp), the weakly fiuctuat- 
ing scaling density 4)2 hi first approximation is associated 
with the density of entropy only, 02 = A(pS). There are 
two important thermodynamic consequences of complete 
scaling that can be checked experimentally. Firstly, the 
"diameter" p^ should contain two non-analytical contri- 
butions, associated with the terms asAP and 62A/}, in 
the scaling fields: 



Pd - 1 = 03 (1 + 0.3) (i>l + b24>2 + 



Di 



AT 



2/3 



D2 



AT 



D. 



AT 



(10) 



where Di = a3PQ/(l + 03) and D2 ~ 62^0/(1 — ") 
with Pq and Aq being the amplitudes in the asymp- 
totic scaling power laws for the liquid/vapor densities. 



Ap = ±Po 



AT 



/3 



An 



AT 



the two-phase region, Cy/fc 
A+/0.523) [5c]. Note, since 2/3 < 1 

2/3 

should dominate near the critical point 



Di 



AT 



and isochoric heat capacity in 

+ ... (Aq-^ 

a, the term 

Sec- 
ondly, the presence of this term implies a so-called Yang- 
Yang anomaly: the divergence of the heat capacity in 
the two-phase region is shared among the second deriva- 
tives of pressure and chemical potential |l2| . Experimen- 
tal verification of complete scaling is a very challenging 
task. The nonanalytical contributions in the "diameter" 
are usually not large enough to be separated unambigu- 
ously. Attempts to fit some experimental and simula- 
tion data to Eq. (10) showed very poor conversions [l3J, 
mainly because of a strong correlation between the linear 

l-Q 



and £'2 



AT 



terms. Experimental tests of the Yang- 
Yang anomaly are even more controversial since traces of 
impurities can easily mimic such an anomaly, thus mak- 
ing any conclusions unreliable |1!tI |. 

We have been able to reliably determine the two asym- 
metry coefficients, 03 and 62, and to conclusively prove 
the validity of complete scaling by combining accurate ex- 
perimental and simulation liquid-vapor coexistence and 
heat-capacity data. We have exploited the fact that the 
coefficients D2 and P3 in Eq. (10) are actually coupled. 
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FIG. 1: Singular diameters in SFb (a) and N2 (b). Exper- 
imental data: SFe [16c] and N2 01 . Curves: solid - fit to 
Eq. (10), dashed - 2/3 term, dotted - 1 — a and linear terms. 
Heat-capacity source [16d] (SFe) and [16i] (N2). 



As the weakly fluctuation scaling density 02 is the crit- 
ical part of the entropy density, in the two-phase region 
at average density p — p^ 






-dT ■ 



1 



AT 



Br 



AT 



(11) 



where Cy is the critical part of the isochoric heat 
capacity and B^ the so-called " critical background" , 
fluctuation-induced analytical part of the heat capacity. 
The critical background can be obtained from a ratio be- 
tween Bci, A^ , and the non- asymptotic heat capacity 
amplitude Af |1^] and by subtracting the "ideal gas" 
from the total heat-capacity background. Both proce- 
dures yield very similar values of i3cr- Since D3 — — 62^cr, 
Eq. (10) contains only two adjustable coefficients, 03 
and 62- We have examined a number of systems, real 
fluids and simulated models [3|, [ig , for which we could 
find both heat-capacity and coexistence data in the range 
AT < 0.01 . In this range the terms of higher-order than 
linear in Eq. (10) are within experimental errors. Exper- 



imental data closer than 



AT 



< 10 were avoided as 



they might be affected by errors in p^ and Tc and by other 
factors, such as gravity, impurities, etc.. For all systems 
studied we have been able to obtain reliable values of 
Bci and conclusively separate two singular contributions 
to the diameter. Two typical examples that represent 
two different kinds of asymmetry in fluid criticality are 
shown in Fig. 1. 

In diameters of some fluids, such as SFg, C2F3CI3, and 



n-CyHie, the 



AT 



2/3 



term dominates (03 is relatively 



large and positive ) while in many other fluids, such as 



FIG. 2: Complete-scaling asymmetry coefficients 03 and 62 
versus reduced critical density p* — p^(8^n). VDW is a 
modified-by-fluctuations van der Waals fluid 19] with a short 
interaction range (R = (p*)^ = 0.5). HCSW is a simulated 
hard core square well model [l^. For a similated restrictive 
primitive model (RPM) [l^] as = 0.14 and 62 = -0.48 (off 
scale) with p* ~ 0.22. The solid curves are given as a guid- 
ance. 



HD, Ne, N2, and CII4, the two singular contributions in 
diameter partially compensate each other (03 is small and 
negative), creating an illusion of rectilinear diameter even 
close to the critical point. In Fig. 2 the two asymmetry 
coefficients are plotted against the dimensional density p* 
defined as p* — P(.(8^o), where ^p is the amplitude of the 
correlation length (representing the range of interactions) 
obtained from the heat-capacity amplitude Aq through 
the two-scale factor of universality, Aq p^^^ = 0.171 [5c]. 
A general trend in the two sources of asymmetry is clear: 

2/3 



the 



AT 



singularity is a dominant contribution into 
the singular diameter if the ratio of critical volume per 
molecule {Pc^) to "interaction volume (S^q)" is large 
0. 

Assuming that the relations between physical fields 
and scaling fields are not affected by fluctuations, we have 
developed a method to obtain the asymmetry coefficients 
03 and 62 from mean-field ("classical") equations of state. 
By combining Eqs. (3) and (7-9), we obtain 



03 



03 



Spii 5^30 fill 5^3Q^;^2 



with p^j — d'^^^ p/dp'dT^ . We have obtained 03 and 62 
for a few classical equations of state, the fine-lattice dis- 
cretization model (crossover between the van der Waals 
fluid and lattice gas) [l8j, the Debye-Hiickel and Flory- 
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We conclude that the asymmetry in near-critical fluids 
originates from two sources: one is a coupling between 
the density and density of entropy, another one is a non- 
linear coupling between the density and molecular vol- 
ume. Both sources can be incorporated into symmetric 
Ising criticality by a proper mixing of physical fields into 
scaling fields. 
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FIG. 3; Crossover diameter in a van der Waals equation of 
state modified by fluctuations witli a short interaction range 
R = (p*)^ = 0.5. Tliick solid curves are the phase boundary. 
(a) two contributions in the singular diameter (solid line): 
1 — a and linear term (dashed line) and 2/3 (dotted line); (b) 
Crossover between rectilinear diameter (dashed-dotted line) 
and singular diameter (solid line) in a broder critical region. 



Huggins models. The coefficient of rectilinear diameter 



D = 



1 






-B'' 



as Mso 






(13) 



where Bq and ACy are mean-field amplitudes of coex- 
istence densities and heat-capacity discontinuity. Close 
to the critical point the rectilinear diameter, affected by 
fiuctuations, splits into two nonanalytical terms, shift- 
ing the critical density In Fig. 3 crossover between 
rectilinear diameter and complete-scaling singular di- 
ameter is shown for the van der Waals equation of 
state renormalized by fluctuations [l9j- A fluctuation 
shift in the van der Waals critical density is mainly 



controlled by the 



AT 



2/3 



singularity since the van der 
Waals value of as is relatively large (as /(I + as) = 
&2/A11 = 0-2) • In the same fashion we have calculated 
crossover between the discontinuity A ( cPfi/dT^ J = 

V / cxc 

(Mn/M3o)[(-M2i/Mii) + (3^40/5^*30)] = -as/(l + 
a3){ACv/k-B), and a divergence known as the Yang- 



a:i[d?P/dT'' 



Yang anomaly, (d?ji/dT' 

— 03/(1 -t- as) A,^ AT . A renormalization-group treat- 
ment of the (Ap)^ term instead predicts a cusp, contain- 

10] . Similarly, 



AT 



-a-p 



mg a term oc 

there should be a term oc as/(l 



AT 



0.865 



as) 



AT 



instead of 



AT [20| in the so-called Tolman's length, a curva- 

ture correction to the surface tension of a liquid droplet. 
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